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Plasmonic resonance of a metallic nanostructure results from coherent motion of its conduction 
electrons driven by incident light. At the resonance, the induced dipole in the nanostructure is 
proportional to the number of the conduction electrons, hence 10 7 times larger than that in an atom. 

The interaction energy between the induced dipole and fluctuating virtual field of the incident light 
can reach a few tenths of an eV. Therefore, the classical electromagnetism dominating the field 
becomes inadequate. We argue that quantum electrodynamics (QED) should be used instead as 
the fundamental theory to describe the virtual field and its interaction with the electrons. Based on 
QED, we derive analytic expressions for the plasmonic resonant frequency, which depends on three 
easily accessible material parameters. The analytic theory reproduces very well the experimental 
data, and can be used for rational design of materials. 
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When light interacts with a metallic nanostructure, its 
conduction electrons may undergo collective oscillations 
driven by the electric field of the light. Known as lo¬ 
calized surface plasmon resonance (LSPR), the collective 
oscillations can be tuned by adjusting the shape, size and 
surrounding medium of the nanostructure, which is at the 
heart of the burgeoning field of plasmonics with potential 
applications ranging from photocatalysis fj], [2j] to optics, 
chemical and biological sensing Q, and photo-thermal 
therapeutics, to name a few 

Although the fundamental physical theory of light- 
matter interaction is quantum electrodynamics (QED) 
sei, QED has not been brought to bear on problems 
in plasmonics. The present research paradigm of plas¬ 
monics is rooted in classical electrodynamics where the 
electromagnetic field of the light is treated classically. It 
is taken for granted that QED correction is too small 
to be relevant in practical plasmonic applications. The 
tremendous success of the classical theory has certainly 
reinforced this notion. However, as we will show in this 
paper that the QED correction cannot be ignored in plas¬ 
monic resonance. In fact, for collective excitations such 
as LSPR, the QED correction could result in an energy 
shift on the order of a few tenths of an eV, well within the 
experimental probes. For nanoparticles, neglecting the 
correction would lead to a size-dependence of resonant 
frequency that is contradictory to experiments. The main 
purpose of the paper is to highlight the importance and 
elucidate the consequence of QED in plasmonics. Focus¬ 
ing on plasmonic resonance of metallic nanostructures, 
we illustrate the origin of the plasmonic energy shift and 
derive analytic expressions for the resonant frequency. 
The theory is compared to available experimental results 
on nano-spheres, nano-rods and nano-plates and shows 
promise as a rapid means for screening materials and 
structures. 

The present theory of plasmonic resonance rests on 


a classical description of the electromagnetic field, ir¬ 
respective of whether the electrons are treated quan¬ 
tum mechanically or not. In QED, the electromagnetic 
field is quantized and there exists a virtual fluctuating 
electromagnetic field in vacuum carrying zero-point en¬ 
ergy. The fluctuating virtual field results in a fluctua¬ 
tion of electronic coordinates, leading to a change in the 
Coulomb energy between the electrons and the positive 
ions. This energy change owing to the virtual electromag¬ 
netic field is known as Lamb shift in atomic physics, and 
usually extremely small 0. In a plasmonic nanostruc¬ 
ture, the wavelength of the visible light, the beam width 
and the skin-depth are typically greater than the size of 
the nanostructure |ll|,[12|, hence all conduction electrons 
in the nanostructure undergo coherent and collective os¬ 
cillations. For a nanoparticle of a radius of 100 ran, the 
coherent oscillation could involve ~ 10 7 conduction elec¬ 
trons, leading to an induced dipole moment that is 6-7 
orders of magnitude stronger than that in an atom. Since 
the interaction energy of the fluctuating virtual field with 
the conduction electrons is proportional to the induced 
dipole, the energy shift can reach several tenths of an eV, 
which can be measured experimentally. 

Let H e be the electronic Hamiltonian of the nanostruc¬ 
ture, and (|a), \b), • • ■ } and {E a , Eb, • • • } as the eigen¬ 
states and eigenvalues of H e , respectively. Let Hf be 
the Hamiltonian of the incident electromagnetic field in¬ 
cluding both the external field and the virtual field. The 
eigenstates of Hf are labeled as |nk,e,, fik 2 e 2 • • •), where 
rtkie! is the occupation number in the photon state kiei 
(ki is the wave-vector and ei is the polarization vector). 
The total Hamiltonian of the system is thus 

H = (H e + H f ) + U, (1) 

where U is the interaction between the conduction elec¬ 
trons and the virtual field. Since the wavelength of the 
field is much larger than the size of the nanostructure, 
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the interaction energy U can be expressed as I,® 

U = -Y, di-E( rj ), (2) 
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where E(rj) is the electric field at position r, y of the jth 
electron, and dj is the dipole operator of the jth electron. 
The extinction spectrum of the nanostructure, comprised 
of the absorption and scattering spectrum of the inci¬ 
dent photon, is the primary physical quantity that can 
be measured experimentally and calculated theoretically. 
According to QED [§], the resonant frequency ui res of the 
extinction spectrum consists of two contributions: 


a^res — Mm T A. 


(3) 


hu> m represents the excitation energy of the plasmons. 
Since the plasmons or the coherent oscillations are the 
eigenstates of H e , denoted by 4' m , is the correspond¬ 

ing eigenfrequency. A is the frequency shift resulted from 
the interaction of the electrons with the virtual field. 

Here we derive an analytical expression for this eigen¬ 
frequency, with the oscillations along one of the major 
axes of the nanostructure. To a good approximation, 
nanostructures such as spheres, rods, and circular plates 
can be modeled as ellipsoids with a uniform volume po¬ 
larization. Suppose that a metallic ellipsoid with a di¬ 
electric function Si = e\ + is" is embedded in a medium 
with a dielectric function £2 = 4 + i £ 2 > and the elec¬ 
tric field of the incident photon, £’® xt , is along the z-axis 
of the ellipsoid. If the longest dimension the ellipsoid, 
tmax, is much smaller than the wavelength A, the beam 
width W of the incident light [ 11 ], and the skin-depth 
S = A(27r£ , 1 /1 ^ 2 ) -1 , the total induced dipole moment in 
the ellipsoid is given by 12. 13j|: 


V lz = a z Ve 0 Et xt , with a z = 


£ 2 (£i - 1 ) 


£2 + (£1 - £ 2 )4 2 ) 


(4) 


Here V is the volume of the ellipsoid; eo is the permeabil¬ 
ity of free space and n^ is the depolarization factor along 
the z-axis. We now make two assumptions: (1) 4 is a 
constant, which is a reasonable assumption for commonly 
used media, such as vacuum, SiC >2 and polyvinyl alco¬ 
hol, in plasmonic applications. (2) The dielectric func¬ 
tions of the nanostructures are primarily determined by 
the conduction electrons 14|, which is also a reasonable 
approximation. Under these assumptions, the dielectric 
functions of the nanostructure can be expressed as: 


where vf = h(3n 2 ri) 1 / 3 /m is the Fermi velocity of the 
metal and L e g- is the effective dimension of the ellipsoid 
along the polarization direction. A is a dimensionless, 
positive constant on the order of one. 70 is the decay 
rate resulting from electron-phonon, electron-impurity, 
and electron-electron interactions [14j . and can be esti¬ 
mated by 70 = ne 2 /mo with a as the conductivity of 
the metal. The second term in Eq.© stems from elec¬ 
tron scattering with the surface, and is the only term 
dependent on the size of the nanostructure. 

In the classical electromagnetism, the extinction cross- 
section rTg X (u;) of the nanostructure is proportional to 

{[4 + (4 - £>W] 2 + [£" + (£" - 4 V *)] 2 }- 1 0. We 

show in the Supplementary Information that the eigen¬ 
frequency M m maximizes <Tg X (wh Therefore M m is the 
root of the following equation [18 h26I ]: 

e 2 + ( £, i ~ e 2) n< ‘ 2 '* = 0- (7) 


Eqn. 0 yields the expression for the eigenfrequency t u m - 




1 + 4 (^- 1 ) 


—) 2 ] 1/2 . 


( 8 ) 


Here and later we drop the superscript z for brevity when 
there is no confusion. 

Based on non-perturbative many-body quantum the¬ 
ory, the frequency shift A can be determined from the 
Green’s function of the Hamiltonian 0. Specifically, 
hA = Re(-R) and R is the shift-operator of the Green’s 
function, given by |9j] 


r |(a; ke|I/|\l/ m ; 0)| 2 

. Em E a Hu) ke 

a ke 


+E E 

ake a' k'e' 


(9) 


(^ m ; 0 |U|a;ke)(a;ke|U|a / ;kV)(a / ;k'e'|U|^ rra ; 0 ) 

(A m E a /zuk e )(-E m E a > licak'e') 

where E m is the energy of the eigenstate |’I' m ; 0) rep¬ 
resents the many-body state at which the incident photon 
is absorbed and the plasmon is excited. |a;ke) denotes 
direct-product state of a many-electron state |a) and a 
virtual photon state with a wave-vector k and a polar¬ 
ization vector e. Finally, we arrive at 


Mm Pz 

4 1 +/V 


( 10 ) 


£1 = 1 - 


SmI 


UJ* 


v 2 ’ 


Sim 2 p 

M 3 + 7 2 M 


(5) 


Here m p = (ne 2 /me o ) 1 ^ 2 is the plasmon frequency, n and 
to denote the electron number density and mass of the 
electron. S(u) represents intra-band oscillator strength 

of the conduction electrons._ 7 is the decay rate of the 

quasiparticles, given by [14|-|l7|: 


7 = 70 + Av F /L eS , ( 6 ) 


with 


v 4 2 
Am M *)] 2 



) 2 }- 


( 11 ) 


Here A m = 27rc/w m - In Eqn. (fill) . e\ and e" are evalu¬ 
ated at m = M m - As /3 Z > 0, the frequency shift A is less 
than 1/4 of M m . 

According to Ean.dSl flOl) . there are three contributions 
to M res - Among them, the dominant one is the first term 
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in Eqn. © since « 1; the second term in Eqn. © 
is the smallest among them. The dominant term depends 
only on the shape of the nanostructure through the depo¬ 
larization factor, rv- z \ Hence the resonant frequency w res 
of the nanostructure is determined primarily by its shape 
as opposed to its size. This fact has been well established 
and exploited in plasmonics ll Hi More importantly, 7 
in Eqn.© is a monotonically decreasing function of L e g 
as indicated in Eqn.©, hence Lo m is a weakly increasing 
function of the particle size. If there were no correction 
term A, the resonant frequency w res (=w m ) would have 
been a monotonically increasing function of the parti¬ 
cle size^ which is opposite to the experimental observa¬ 
tions Bii [25l l . The QED correction term, A, as 
a decreasing function of the particle size, reverses the in¬ 
correct size-dependence of the classical electromagnetic 
theory and renders w res consistent with the experiments. 
The failure of the classical theory has also been discussed 
by Scholl et al. 21] who attributed the opposite size- 
dependence to the inappropriate use of macroscopic di¬ 
electric functions in the nanoparticles. The macroscopic 
dielectric functions failed to capture the effects of discrete 
energy levels and the fact that only certain electronic or 
plasmonic transitions are allowed in the nanoparticles. 
To remedy the classical theory, Scholl et al. proposed 
a phenomenological model based on discrete energy lev¬ 
els. Although the model yielded an improved agreement 
to the experimental data, it did not consider the quan¬ 
tum effect of the electromagnetic field. As a result, the 
model cannot guarantee the monotonically decreasing 
size-dependence of w res , as revealed in experiments and 
the present theory. Nonetheless, Scholl’s model is valu¬ 
able contribution and could be combined with the present 
theory to form a more comprehensive microscopic picture 
of plasmonic resonance. 


If the electric field of the incident light is perpendicu¬ 
lar to the z-axis, the depolarization factor in the normal 
direction has to be worked out. We have derived the 
corresponding equations, which are included in the Sup¬ 
plementary Information. Moreover, if the electric field 
of the incident light is along an arbitrary direction, the 
total induced dipole is a vector sum of the components in 
each major axis iB Finally, for a spheroid with its 
rotational axis along z, the depolarization factors n^ x \ 
n (y), n (z) anc i the resonant frequency w| es (E ext ||z) and 
(E ex t -Lz) can be calculated analytically as well. For 
a general ellipsoid, the corresponding quantities have to 
be evaluated numerically. 

To validate the proposed theory, we apply it to various 
metallic nanostructures including nano-spheres, nano¬ 
rods, and nano-plates. First, we examine the size- 
dependence of w res in nano-spheres. Since all spheres 
have the same shape or the depolarization factors (n = 
1 /3) [ijJ, the nano-spheres of the same metal would yield 
the same w m for a given surrounding medium. Accord¬ 


ing to Eqn. m El A depends only on the volume of 
a sphere, thus w res is a monotonically decreasing func¬ 
tion of the sphere diameter D. In Fig. 1(a), we compare 
the theoretical prediction to the experimental data taken 
from Fig. 3(b) of reference [21] for Ag nano-spheres. The 
dielectric constant of the surrounding medium e' 2 is 1.69 
as measured in the experiment [2lj| . We find that the 
theoretical prediction agrees very well to the experimen¬ 
tal data as long as the two fitting parameters A and S 
are chosen reasonably, in this case A = 0.03, 5=1. As 
displayed in Fig. 1(b), the size dependence of w res is en¬ 
tirely contained in A while cu m is essentially flat. Thus 
the size-dependence of the nano-spheres originates exclu¬ 
sively from the quantum nature of the electromagnetic 
field! 


Second, we compare the theoretical prediction to the 
experimental results [TD for Au nano-rods embedded in 
silica (e' 2 = 2.15). In Fig. 2 (a), the experimental res¬ 
onance frequency cu| es as function of the length L and 
diameter D of the nano-rods is shown in circles, while 
the theoretical prediction is the surface. In this case, the 
two fitting parameters are A = 0.6, S = 1.4. Because 
L ~ 32-70 nm [19], larger than the size of the nano¬ 
spheres, the electron scattering at the surface becomes 
more important, thus A is larger. For the similar reason, 
the oscillator strength S is also larger than the nano¬ 
spheres, as discussed in the Supplementary Information. 
There is an overall good agreement between the theory 
and experiment, down to the size of 8.5 nm [19]. The 
theoretical prediction for the polarization direction per¬ 
pendicular to the nano-rod axis is displayed in Fig. 2 
(b). 


We next demonstrate the validity of u > m , which cannot 
be measured directly by experiments. Hence, we compare 
the theoretical prediction of u) m to a set of computational 
results obtained from time-dependent orbital-free den¬ 
sity functional theory (TD-OFDFT) simulations 271 for 
a Na nano-rod embedded in vacuum (£ 2 = 1 ). A number 
of (L,D) combinations including (5.79,0.86), (5.79,1.41), 
(5.79,1.93), (5.79,2.23), (5.79,3.54), (5.79,4.76), and 
(5.79,5.46), in the unit of nm, are considered. The two 
fitting parameters are A = 0.6, S = 0.8 for E || axis 
and A = 0.6, S = 0.92 for E _L axis. There is an excel¬ 
lent agreement between the theoretical predictions and 
the TD-OFDFT results for both polarization directions 
as shown in Fig. 3, which validates the derivation of w m . 
In Fig. 2 and 3, one may notice that Huires(L , D ) has 
an opposite size-dependence as hu>^ s (L, D ), owing to the 
opposite (L,D) dependence of n^ and n^ x \ 


Third, we switch to a Au nano-plate whose rotational 
symmetric axis is along z. The theoretical predictions for 
uj^, s (D) are compared to the experimental results fl 8 [ . 
The circular nano-plate has a height L = 20nm embed¬ 
ded in a medium with a refractive index of 1.26. The 
fitting parameters A = 0.6, S = 0.86 yield an excellent 
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agreement between the theoretical predictions and the 
experimental data as shown in Fig. 4. Similar agree¬ 
ment between the theory and experiments is also ob¬ 
served for A1 and Pt nano-plates which is presented in the 
Supplementary Information. The opposite D-dependence 
between Wres(-D) and to^ s (D) is due to the opposite D- 
dependence of their respective depolarization factors. 

QED could have more profound implications in plas- 
monics than what is presented in this paper. For ex¬ 
ample, it is known that when a nano-antenna is placed 
next to a metallic nanostructure, there is an interaction 
between the nano-antenna (an emitter) and the virtual 
field. The interaction could c hang e the directional ra¬ 
diation pattern of the antenna 28]. analogous to cavity 
QED [29|. The present work, on the other hand, focuses 
on the interaction between an absorber (the plasmonic 
nanostructure) and the virtual field. Such interaction 
could also change the induced magnetic moment of the 
metallic nanostructure, as well as the polarization of the 
incident light. 

To summarize, we propose that QED is crucial to un¬ 
derstand the plasmonic resonance in metallic nanostruc¬ 
tures. The coherent motion of the conduction electrons 
in the nanostructure could lead to a large induced dipole 
moment, which interacts with the virtual field and re¬ 
sults in a significant shift in the resonant frequency. The 
frequency shift is the key to reconciling the theoreti¬ 
cal predictions and experimental observations on size- 
dependent plasmonic resonance. Based on QED, we have 
derived analytic expressions for the plasmonic resonant 
frequency, which depends on three easily accessible ma¬ 
terial parameters - the dielectric constant e' 2 of the sur¬ 
rounding medium, the number density n of electrons and 
the conductivity a of the metal. The analytic expres¬ 
sion are shown to reproduce very well the experimental 
data for nano-spheres, nano-rods and nano-plates, and 
can be used readily for estimating the resonant frequency 
of plasmonic nanostructures as a function of their geom¬ 
etry, composition and surrounding medium. 

Supplementary Information is linked to the online 
version of the paper. 
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FIG. 1: (a) The plasmonic resonance energy &j res as a 

function of the nano-sphere diameter D, compared between 
the theory (solid curve) and experimental data (triangles). 
(b)The contribution of w m (dot) and |A| (dash) as a function 
of D. 
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FIG. 2: The plasmonic resonance energy of a Au nano-rod 
embedded in Si02 as a function of its diameter D and length 
L, compared between the theory (surface) and experiment 
(circle). The polarization direction is parallel (a) and perpen¬ 
dicular (b) to the nano-rod long-axis. 
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FIG. 3: The eigenfrequency w m as a function of L/D for a Na 
nano-rod in vacuum, determined from the theory (curve) and 
the TD-OFDFT calculations (triangles). The polarization di¬ 
rection is parallel (a) and perpendicular (b) to the nano-rod 
long-axis. 
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FIG. 4: The plasmonic resonance energy as a function of di¬ 
ameter D for a circular gold nano-plate with a fixed height L 
= 20 nm. The results from theory are shown in solid curves 
while the experimental data is shown as triangles; The polar¬ 
ization direction is parallel (a) and perpendicular (b) to the 
long-axis of the plate. 














